Estimation of Persistence Lengths of Semiflexible Polymers: Insight from Simulations 
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The persistence length of macromolecules is one of their basic characteristics, describing their 
intrinsic local stiffness. However, it is difficult to extract this length from physical properties of the 
polymers, different recipes may give answers that disagree with each other. Monte Carlo simulations 
are used to elucidate this problem, giving a comparative discussion of two lattice models, the self- 
avoiding walk model extended by a bond bending energy, and bottle-brush polymers described 
by the bond fluctuation model. The conditions are discussed under which a description of such 
macromolecules by Kratky-Porod worm-like chains holds, and the question to what extent the 
persistence length depends on external conditions (such as solvent quality) is considered. The 
scattering function of semiflexible polymers is discussed in detail, a comparison to various analytic 
treatments is given, and an outlook to experimental work is presented. 
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I. INTRODUCTION 

Flexibility of chain molecules (or lack of flexibility, 
respectively) is one of their most basic general proper- 
ties [H-Q . It affects the use of macromolecules as building 
entities of soft materials, and controls some aspects of the 
functions of biopolymers in a biological context. Thus, 
it is important to understand its origin in terms of the 
macromolecular chemical architecture, and the extent to 
which it depends on external conditions (temperature, 
solvent quality if the polymer is in solution, as well as 
polymer concentration), and one therefore needs to be 
able to characterize macromolecular flexibility or stiff- 
ness precisely. The quantity that is supposed to describe 
the local intrinsic stiffness of a polymer is termed "per- 
sistence length" and often it is introduced (e.g. d, Q) 
as a length describing the exponential decay of orien- 
tational correlations of segments with the length of the 
piece of the chain separating them. Thus, let us con- 
sider a linear macromolecule composed of segments vec- 
tors {Hi, i = 1, - ■ ■ , N)}, all having the same bond length 

((af) — ^1, if we wish to allow for thermal fluctuations 
of the length of these segments). Then it is assumed that 
the correlation of two segments i, j , that are s = \i — j\ 
steps along the chain apart, varies as 

(cos0(s)) = (a; • aj)/{af) — exp{~sib/ip) , s — ?► oo, 

(1) 

where £p is the persistence length. 

In fact, Eq. ^ holds for models of linear polymer 
chains that strictly follow Gaussian statistics (for large 
distances between monomeric units), however, Eq. ([l} is 
not true for real polymers, irrespective of the considered 
conditions: for dilute solutions and good solvent condi- 
tions one rather finds a power law behavior Q 

(cos 6l(s)) cx fi-^ , P^2{l^iy), 1< s < TV . (2) 

Here v is the well-know Flory exponent, describing the 



N 

scaling of the end-to-end distance R — J2 '^i with the 

1=1 

number N of segments, (i?^) cx A^^", with ly ss 3/5 
(more precisely [3], = 0.588) in d = 3 dimensions [ll-Q. 
Polymer chains in dense melts do show a scaling of the 
end-to-end distance as predicted by Gaussian statistics, 
(i?^) oc N (i.e., I' takes the mean-field value vmf = 1/2), 
and hence it was widely believed, that Eq. ([T]) is useful 
for polymer chains under melt conditions. However, re- 
cent analytical and numerical work ^ |9i] has shown that 
this assertion is completely wrong, and there also holds 
a power law decay, though with a different exponent, 



(cos e{s)) oc s"'''/^, K s < iV. 



(3) 



More recently, it was also found by approximate ana- 
lytical arguments [loj . and verified in extensive simu- 
lations [ll| that Eq. (jSj also holds for chains in dilute 
solutions at the Theta point. In practice, since asymp- 
totic power laws such as Eqs. ^ hold only in the 
intermediate regime 1 <C s ^ TV and hence one must 
consider the limit TV — > oo, one easily could be misled if 
data for (cos0(s)) are considered for insufficiently long 
chains. As an example Fig. [T] presents simulation re- 
sults for the simple self-avoiding walk (SAW) model on 
the simple cubic (sc) lattice, where an attractive energy 
e between neighboring occupied sites (representing the 
effective monomers of the chain) occurs and the temper- 
ature is chosen as ksT /e = 3.717 which is known to re- 
produce Theta point conditions for this model jl^ . One 
can see clearly that the data for — oo and s > 10 do 
approach Eq. but for finite N systematic deviations 
from Eq. ([3]) clearly are visible already for s = N /lO. 
On the semi- log plot, for rather short chains one might 
be tempted to apply a fit of an exponential decay pro- 
portional to exp(— s£fc/£p) to the data for rather large s, 
but resulting estimates for ip/ lb are not meaningful at 
all: for the considered model, the chain is fully flexible, 
any reasonable estimate for tp/lb that describes the local 
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FIG. 1. Semi-log plot (a) and Log-log plot (b) of (cos^(s)) 
versus s as obtained from Monte Carlo simulations (as de- 
scribed in [lU l using the pruned-enriched Rosenbluth method 
(PERM algorithm for a self-avoiding walk with nearest- 
neighbor attraction e, under Theta point conditions. The full 
curve in (a) and straight line in (b) represents the relation 
(cos6l(s)) = 0.16s-^''^ 



intrinsic stiffness of the chain should be (i) of order unity 
(see Fig. [l}i, £p/£b ~ 0.94), and (ii) independent of N. 
Both conditions are dramatically violated, of course, if 
estimates for £p/£b were extracted from fits to an expo- 
nential decay in this way. 

Since the intrinsic stiffness of a chain is a local prop- 
erty of a macromolecule, one might alternatively try the 
recipe to either fit Eq. ([T]) in the regime of small s to the 
data, or assume that Eq. ([Ij holds for s = 1 already and 
hence 

V4 = -l/ln((cos0(l))) . (4) 

This recipe works in simple cases, such as the SAW 
model where an energy et associated with bond bend- 
ing is added (every kink of the walk by ± 90° on the 
sc lattice costs £&), see Fig. [2l but it fails for molecules 
with more complex chemical architecture, such as bottle- 
brush molecules [11-13. The dramatic failure of Eq. 
for bottle-brush polymers is understood in terms of their 



FIG. 2. Semi-log of (cosS(s)) vs. s for a semifiexible version 
of a SAW model on the sc lattice (a), cf. text, and the bond- 
fluctuation model of bottle-brush polymers under very good 
solvent conditions [O, [l3l (b). Part (a) refers to the chains of 
length N = 50000, and several choices of the parameter qt, = 
exp{—eb/kBT) controlling the chain stiffness, namely qt = 
0.4, 0.2,0.1,0.05 and 0.02. Using Eq. @, the straight lines 
indicate the exponential decay exp{—s£b/£p) for the choices 
of qb- Part (b) refers to the case of bottle-brush polymers 
where every effective monomer of the backbone has one side 
chain of length Ns — 24 grafted to it, and several choices of 
backbone chain length Nb- Here £p/£b = — 1/ ln((cos 6(1))) 
has been extracted from the chain backbone only. 



multiscale structure (Fig. ^ : The side chains lead to a 
stiffness of the backbone on a mesoscopic scale, even if on 
the local scale of nearest-neighbor bonds the backbone 
is still rather flexible. The question of understanding 
this stiffening of bottle-brus h polym ers because of their 
grafted linear side chains [U ITSmSj l or grafted branched 
objects [sgl - li^ is an issue of longstanding debate in the 
literature. 

Complex polymer architecture is only one out of many 
reasons which make the analysis of bond orientational 
correlations based on Eqs. ([l]) or (|4l) problematic. In 
dilute solutions we expect that a nontrivial crossover oc- 
curs when the solvent quality is marginal, i.e. close to 
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FIG. 3. Sketch of the muhiple length scales that one may 
define for bottle-brush polymers (schematic): While on scales 
where one resolves the effective monomers of both the back- 
bone and the side chains, the correlations of backbone vec- 
tors a' (which have length £(,) and side chain bond vectors af 
(which have length t^s) can be studied, as well as end-to-end 
distances R^.a of side chains and of the backbone Re,bb (and 
the corresponding gyration radii). In the bond- fluctuation 
model of Fig. [2] for simplicity no chemical difference between 
backbone and monomer was considered, so It — Is was cho- 
sen. The "microscopic" contour length of the backbone then 
is Lh = Nblb, if the backbone has A^i, bonds. On a coarse- 
grained level the bottle-brush resembles a worm-like chain 
of thickness (cross-sectional radius) Res and contour length 
Lcc < Lt, which is locally straight on the scale of the persis- 
tence length Ip. 



the Theta point a large size of "thermal blobs" [i^ 
exists, such that for values of s along the backbone of 
the chain corresponding to distances r(s) > one ex- 
pects that excluded volume effects are visible and hence 
Eq. ([2]) should hold. For semidilute solutions 43] , on the 
other hand, in the good solvent regime the inverse ef- 
fect occurs: there exists a screening length ^(c) depend- 
ing on the polymer concentration c (also called size of 
"concentration blobs" such that excluded volume 

effects are pronounced for r(s) < ^(c) but are absent 
for r(s) > ^(c). Then Eq. ^ holds for the latter case 
and Eq. ([2]) for the former, for rather flexible chains. If 
the chains are semiflexible, in favorable cases (e.g., for 
simple chemical architecture of the polymers) we might 
observe Eq. ([T]) for I < s < s* where s* depends on the 
local intrinsic stiffness of the chain, which we wish to 
characterize by ip. Then the question arises whether s* 
is smaller than any of the other crossover chemical dis- 
tances (due to marginal solvent quality, described by a 
Flory-Huggins parameter x with (1/2 — x) ^ f [il,E3|j 
or due to nonzero c) or not. The conclusion of this dis- 
cussion is that the behavior of bond orientational correla- 
tions {(ii ■ dj) is subtle, and not always suitable to obtain 



straightforwardly information on the intrinsic stiffness of 
macromolecules; as a further caveat we mention that in 
general it is also not true that this correlation depends 
on the relative distance s — \i — j\ only: it matters also, 
if one of the sites is close to a chain end. 

Another popular definition is the local persistence 
length £p{i) defined as [H, 0] 

£p(z)/4 = (a,.i?)/(a2). (5) 

However, it has been shown by renormalization group 
methods that in good solvents one has, for N oo, 
£p{i) cx [i{N - i)]'^"^^, i > f, so the behavior of £p{{) 
in the chain interior clearly is unsuitable to conclude 
anything about the local stiffness of a chain under good 
solvent conditions, and this conclusion has been corrob- 
orated by simulations [Tf|, Sometimes it has been 
argued that a better choice is to take the correlation 
between the first bond vector and the end-to-end dis- 
tance, ip{l) [46j . However, since in a macromolecule the 
chemical nature of the end monomer always differs from 
inner monomers, one can never expect that ip{l) pre- 
cisely characterizes the local stiffness of a linear macro- 
molecule in the inner parts of a chain. Moreover, since 
(i?^) reflects all the crossovers (due to "thermal blobs" 
etc.), [i^, as discussed above, it is premature to expect 
that ip{l) stays unaffected from them. We also note that 
for d = 2 dimensions under good solvent conditions it has 
been shown [d^l that £p{l) oc InTV — )■ oo as iV ^ oo, so 
in this case £p{l) clearly is not a useful measure of the in- 
trinsic stiffness of a chain at all. Since Eq. ([5]) is difficult 
to extract from any experiments, and inconvenient for 
simulation studies due to high sampling effort, we shall 
not discuss Eq. ^ further in the present paper. 

Experimental studies try to extract the persistence 
length either from scattering analyses of the single chain 
structure factor (e.g. [28l. [33. l33438i , idS] ) or from analy- 
ses of extension versus force measurements of stretched 
chains (e.g. [49l - [57| ). However, the interpretation of the 
latter experiments must rely on a theoretical model of 
the extension versus force curve. While this task is sim- 
ple for ideal random walk models of polymers [3, d, [s^ 
and also for semiflexible polymers when excluded vol- 
ume is neg lected fsoj . so that the Kratky-Porod (K-P) 
model |6C| of worm-like chains can be used, it is very 
difficult (due to multiple crossovers [6l|,[63]) if excluded 
volume effects are included. These excluded volume ef- 
fects cause an intermediate nonlinear variation of the ex- 
tension versus force curve (the chain is then a string of 
"Pincus blobs" ^635), making the estimation of the per- 
sistence length difficult [62], and this behavior has also 
been verified in recent experiments [56l . [57| . Since we 
have given a recent extensive discussion of this prob- 
lem elsewhere [62], we shall not dwell on this problem 
here further, and focus on the problem how the persis- 
tence length shows up in the single chain structure fac- 
tor S{q). Here the key idea is that the scattering in- 
tensity S{q) at scattering wavenumber q yields informa- 
tion on the structure of the macromolecule at a length 
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scale A = ^ir/q. This problem also is subtle, even in the 
framework of simple models (see Fig. |3^,b,c,d) used for 
simulations. If A is of the scale of the cross sectional 
radius Res for the models (a,b,d) or the lattice spacing 
in (c), local structure on the scale of effective subunits 
is revealed: soft (a) versus hard (b) effective cylinders, 
hard spheres in (d), but one could also conceive a chain 
where soft spheres are jointed, etc. When one considers 
semiflexible chains with no excluded volume, the persis- 
tence length £p would be just one half of the step length 
in cases (b), (d), where one then requires a strong 
bond angle potential to make these chains semiflexible 
rather than flexible; however, as emphasized above, such 
models neglecting excluded volume completely will in- 
evitably imply Eq. ([T|), which is inappropriate for real 
polymers under all physically possible conditions. So the 
information on chain stiffness, as described by the per- 
sistence length, is hidden in some intermediate range of 
wavenumbers. E.g., for the model (c), which will be used 
extensively in the rest of the paper (but in d = 3 di- 
mensions, since the case of d = 2 is rather special 62 1 
as will be discussed below), we need wavenumbers in the 

range 27r/ \J {Rf) ^ <? ^ 27r/a, where a is the lattice 

spacing. The aim of the present paper is to present a 
discussion of how one can obtain detailed information 
on intrinsic chain stiffness from the gyration radius of 
the macromolecules and from the structure factor S{q) 
in the suitable intermediate range of wavenumbers q. 

The outline of our paper is as follows: in the next 
section, we summarize some pertinent theoretical results 
on S{q). In the third section, our Monte Carlo simula- 
tion methods are briefly described. In the fourth section, 
a comparative discussion of simulation results for two 
models is given, the bond fluctuation model of bottle- 
brushes (c.f. Fig. |4^), and the self-avoiding walk model 
on the simple cubic lattice with variable bending energy 
(of. Fig. I3J;). The final section contains our conclusions. 



II. SOME THEORETICAL RESULTS ON THE 
STRUCTURE FACTOR OF ISOLATED 
MACROMOLECULES IN SOLUTION 

We consider here a single macromolecule with linear 
chain architecture, assuming a sequence of iV -f- 1 (effec- 



tive) monomeric units at positions r}, j = 1,2,- 
with effective bond vectors aw = 
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1 



■■■ ,N. 



We have in mind application to standard polymers like 
polystyrene (disregarding here the scattering from the 
side groups that are attached to the backbone of the 
chain, see e.g. Rawiso et al. [i^ for a discussion of this 
problem in an experimental context). We also have in 
mind application to the scattering from the backbone 
of bottle-brush polymers (this is experimentally directly 
accessible from neutron scattering '28'] if selective deuter- 
ation only of the backbone is used, while in the case of 
deuteration of the whole macromolecules [s^, 1131 tbis in- 
formation can be inferred only indirectly). Due to the 
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h =2RC! 



FIG. 4. Various models of semiflexible polymers, as discussed 
in the context of simulations. Case (a) shows the snapshot 
picture of a typical conformation of a simulated bottle-brush 
polymer using a backbone chain length Ni, — 1027, side chain 
length Ns = 24, projected into the a;t/-plane (this model is 
discussed in more detail in Sec. 3). Case (b) shows a model 
of freely jointed cylindrical rods of Kuhn step length £k and 
diameter D — 2Rcs, with Res the cross-sectional radius (if 
Res = this leads to a simple off-lattice random walk con- 
figuration, while excluded volume interaction is introduced if 
overlap of the cylinders is forbidden). Case (c) shows the SAW 
model on the square lattice with lattice spacing a (_D = a in 
this case), where 90° bends cost an energy Sb ^ fcsT, so 
the chain consists of straight pieces where ristr steps go in the 
same lattice direction, with Ustr ^ 1. Case (d) shows a model 
of tangent hard spheres with radius Res (and £k ~ 2Rcs)- 



restriction to "effective monomeric units" rather than 
talking about the scattering from individual atoms with 
the appropriate scattering lengths, we clearly disregard 
information on the scale of the length of an effective 
bond, but we then need not discuss experimental prob- 
lems such as contrast factors between the scattering from 
the macromolecule and the solvent ■ The effect of the 
cross-sectional structure of the chain (finite chain thick- 
ness D) is not explicitly considered as well (experimen- 
tally this problem often is approximated in terms of the 
Guinier [73| approximation, writing the observed scat- 
tering intensity Sohs{q) = S{q) exp{—q^Rl/2), with Rc 
some "effective" cross-sectional radius of the chain pisj . 
Thus only wavenumbers qD <C 27r are physically mean- 
ingful: in the case of the lattice model. Fig. HJ:, D — a, 
of course. The structure factor then is defined as 



{N + iy 



N + l N+1 



1)2 (II ^exp[iq-( 



j=l fc=l 



(^1 - rk 



), (6) 



and does not depend on the direction of the scattering 
wavevector q. In d = 3 dimensions, it has the small q 
expansion 

^(g) = 1 - (i?2)gV3 + • • • , q^O, (7) 
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where the mean square gyration radius {R"^) enters 

Af+l Af+1 

(^1) = WTW{Y. Eirl-nf). (8) 



j=i k=j+i 



Other characteristic lengths of the chain molecule are the 
mean square end-to-end distance 



<«') =((!:".)) 

1=1 

and the contour length 



(9) 



(10) 



but neither of these lengths can be inferred directly from 
the scattering. For chains in dense melts or in dilute 
solutions under Theta conditions one typically uses an 
ideal chain approximation (disregarding, e.g., logarithmic 
corrections at the Theta point 3,, J^, ^64j] ) 



(11) 



with Coo a characteristic constant llM5|. In this case one 
introduces an equivalent freely jointed Kuhn chain with 
the same contour length, (i?^) = nki\, where uk is the 
number of equivalent Kuhn segments and Ik their length. 



= Cootb , riK = N/C^ 



N oo. 



(12) 



For a semiflexible worm-like chain with Coo ^ 1 Eq- (1) 
holds. However, since under Theta conditions (and 
melts) Eq. (jl2p is approximately true, one finds 



ip = i{Rl)/{Nlb) , (R^) = 2Vfc^ , N 



oo . 



(13) 



if the relation £p = £k/2 = Coo^6/2 then simply is taken 
as an alternative definition of a persistence length. For 
the simple SAW model of Fig. [T] this gives ip = 0.94 lat- 
tice spacings: but as expected, using this value in the 
simple exponential exp{—s£b/ip) one does not obtain a 
description of the actual data in Fig.[l]on the basis of this 
description, because the actual behavior of bond orien- 
tational correlations is a power law decay, Eq. ([3]). Note, 
however, that the relation £p = ^if/2 makes only sense 
for semiflexible chains for which Coo ^ 1 at the Theta 
point, which is not the case for the model of Fig. [TJ 

In the case of good solvent conditions excluded volume 
interactions invalidate Eq. ([TT|) and one finds instead [1, 

asm 



{Ri 



N ■ 



_ (14) 
with [3| ~ 0.588 instead of the mean field value umf = 
1/2 that appears in Eq. Note that we have defined 

the prefactors of the relations (i?^) cx N"^" , (i?^) ^ ]Sf2-^ 
in Eq. in complete analogy with Eq. 11], but 

we shall see shortly that the lengths £p, £p^ do not play 



the role of a persistence length that describes the local 
intrinsic stiffness of the chains. 

For a better understanding of this problem, in partic- 
ular when £p is very large, it is of interest to consider the 
crossover from the rod limit (that occurs for L < £k, i.e. 
nx < 1) to the Gaussian coil limit. This problem can be 
worked out easily for various models of discrete chains [3- 
Q as well as for the Kratky-Porod model. Describing the 
chain by a continuous curve r(s), s being the curvilinear 
coordinate along the chain contour, the potential energy 
of a particular conformation of the chain is given by 



n = 



f/( 



ds, kBT£p {d^'d). (15) 



In Eq. (IT5|) it is clearly assumed that k is a constant, 
independent of the contour length L (or chain length N , 
respectively), and the same holds for £p. The physical 
interpretation of k is in terms of the local bending stiff- 
ness of the chain. Formula (IT5|) can be used for arbitrary 
values of the ratio L/£p = Up, and one can show i65[ 



{R^ 
2£pL 



1 — exp(— Tip) 



and 



3(i?; 



oL 



3 6 6 r 
1 + — - — 



1 — exp(— Tip) 



(16) 



(17) 



One immediately recognizes that for Hp — 2nK — > oo one 
recovers Eq. ([T3|) . while in the opposite limit the results 
for rigid rods of length L are obtained. 



= 12(i?2) rip < 1 



(18) 



However, the generalization of these results to the good 
solvent case, where excluded volume matters, is not 
straightforward. Of course, for np <C 1 excluded volume 
is irrelevant, Eq. (ITSt remains valid. It turns out, how- 
ever, that in d = 2 Eqs. p^ . P7)) are not valid at all, one 
has no regime of Gaussian chain behavior as described 
in Eq. (1131) . and rather near np = 1 a crossover from 
rigid rod behavior to the behavior of two-dimensional 
self- avoiding walks occurs ^62, i66| (z/ = 3/4) 



(19) 



For d = 3, however, Eqs. (fT6|) . (ITTI) for semiflexible chains 
remain valid for Up < n*{£p) where n*{£p — >■ oo) — )■ oo. 
This crossover contour length L* = n *£p has first been 
estimated by a Flory argument as [gl, IgtI . l68j 



L* oc £l/D^ , n; oc {£p/Df 



(20) 



Note, however, that Flory arguments imply v = 3/5 in 
d = 3 (rather than the precise value v w 0.588 [7]) and 
cannot predict any prefactors in Eq. (I^Ul) ; they are based 
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on a crude balancing of the elastic energy of chain stretch- 
ing (taken as Gaussian) and a mean field estimate of bi- 
nary interactions: having in mind a model description as 
in Fig. |3}d, one takes the second virial coefficient propor- 
tional to the rod volume on the scale of the persistence 
length, V2 (X £pD, and in this way the effective chain di- 
ameter D enters the estimated Eq. (1^0]) (6^ WA [68| . Nu- 



merical results, however, seem to suggest that rather [62 1 
n* cx (ip/D)'' with an exponent C ~ 1.5. 

In any case, the conclusion of this discussion is that 
for semiflexible chains in d = 3 the mean square radii 
as a function of the reduced contour length Up = L/ip 
exhibit two successive crossovers, from rods to Gaussian 
coils near Up — 1 and from Gaussian coils to swollen 
chains (described by Eq. ([T4l) ') near Up = n*. These 
two crossovers have in fact been seen nicely in both ex- 
perime nt 16911 and computer simulation of the model of 

Fig. St !6il7o,[n|. 

We now turn to a discussion how these behaviors show 
up in the scattering function S{q) at larger wavenumbers, 
when Eq. ([7]) does not hold. In the regime rip < 1, when 
the chain behaves like a rigid rod, one can work out the 
scattering function in the continuum limit as (72j 



qL 



qL 



dx 



sin a; 1 — cos{qL) 



qL 



(21) 



while for a discrete chain of + 1 scatterers along a rod 
of length L = N£i, one has 



Srodiq) = 



1 



N 



N + 1 



- 1 



N 



fc=0 



qhk 



< 27r. 



(22) 



It is noteworthy to recall that the large g-limit of Eq. (PT|) 
contains information on the contour length L and shows 
a 1/g decay. 



Srodiq ^ oo) = Tr/iqL) 



(23) 



In the Gaussian regime, that applies for chain lengths 
that correspond to 1 <C np < np{£p) in 3, the struc- 
ture factor S{q) is described by the well-known Debye 
function, 



^Dobyc(g) = |{l - ^ [l - CXP(-X)] } , q^Rl) . 

(24) 

For smaU X, Eq. ([Ml) reduces to Eq. ([7]), as it must be, 
while for large X Eq. (gH) yields S'Dobyc(g) ~ 2/X = 
2/[g^(i??)]. While for flexible chains at the Theta point 
Eq. ([24| is expected to hold for large q, up to qh of 
order unity where effects due to the local structure of 
monomeric units comes into play, the validity of Eq. (j24p 
for semiflexible chains is much more restricted, since then 
the rod to coil crossover matters also with respect to the 
intrinsic structure of these polymers, as it is probed by 



S{q). In oder to discuss this problem, it is useful to 
cast S{q) in the representation of the so-called Kratky 
plot [73], qLS{q) is plotted as a function of qL — Y. 
For rigid rods, one simply would have a linear increase 
of qLS{q) with Y for small F, which smoothly crosses 
over near F = 1 to a flat plateau (which has the value 
TT, cf. Eq. p3p ). For chains where intermonomer dis- 
tances follow Gaussian distributions, at all scales, the 
Kratky plot exhibits a maximum at i^max, and then a 
crossover to a decay proportional to Y~^ occurs. To lo- 
cate this maximum, it is convenient to write qLS{q) as 

\fX(Ll ^ {R'^)) Suchyciq) as a function of X, noting that 
the maximum occurs at X^ax ~ 2.13, i.e. the Kratky 

and 



plot has its maximum at Y,j, 



V2J3L/ 



the height of this maximum also is of order L/y (i?^). 

Using now Eq. (fT3|) in the form {Rg) = L£p/3, we recog- 
nize that the maximum of the Kratky plot occurs at 



iqL)r 



^^{Lltpf^ , 1 < Ljt 



p < np{ip) 



(25) 



and also the height of this maximum scales proportional 
to y/L/ip. However, while for flexible chains under Theta 
conditions (for which £p and £b are of the same order), 
one observes on the Kratky plot for F >> Fmax a decay 
qLS{q) oc q~^, for semiflexible chains a crossover from 
this decay to the plateau value tt (given by Eq. (PSI) ) is 
expected. This is also true for semiflexible chains under 
good solvent conditions, if the persistence length £p is 
large enough so that np{= L/ip) < np{£p), and hence ex- 
cluded volume effects still can be ignored. The descrip- 
tion of this decay of the structure factor from its peak 
towards this so-called "Holtzer plateau" ^74^ has been a 
longstanding problem in the literature (75i - i96i | . Only in 
the limit N ^ oo a simple explicit result derived from 
Eq. ^ is available [S^, 

LqS{q) = TT + -[^qip^ , L oo , q^ oo , (26) 

but we should keep in mind that the limit g — > oo is well- 
deflned for a simple mathematical continuum model such 
as Eq. (fT5|) , while for real chains (and for simulations) the 
regime q£b > 2tt is not at all meaningful. Although the 
decay S{q) oc q^^, that Eq. (1^ predicts for q£p ^ 1, 
is compatible with the power law decay of the Debye 
function, Eq. ([M)) . for large q, 



S{q) 



q^Rl) 



qLS{q) « 6{q£p)- 



g^oo, (27) 



the prefactor in Eq. ([271) by a factor of 9 larger than the 
prefactor of the q~^ term in Eq. ([26]) . so both Eqs. ([24]) 
and (|26p are inconsistent with each other. This inconsis- 
tency is due to the fact that Eq. ([26)) is only accurate for 
q£p > 3, it should not be used for sma ll q£r, ■ After many 
less successful attempts, Kholodenko [86l - l89| achieved a 
description which interpolates between the limiting cases 
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of rigid rods and of Gaussian coils, capturing the scat- 
tering law of both limits exactly, but deviating from the 
exact result ( "exact" refers to the Hamiltonian Eq. ([T5)) , 
so no excluded volume effects are being accounted for) 
in the intermediate regime; this exact behavior is known 
from systematic expansions j93 - l96j whose use requires 
heavy numerical work, and will not be considered here. 
Recently we have shown [o^] that the exact method of 
Stepanow [H, HI] deviates only very little from the ap- 
proximation of Kholodenko \89\ , which can be cast in the 
form 



Siq) = - [/i(a;) - -hix) , x = 3X/2£p , 



(28) 



where 



Ux) = J dzz^-'fiz), 



f 1 sinh(£z) „ < Q/O/- 
f(^\^)E sinhz ' y^i-J/^-^pj 

[ E' sinhz ' y>'J/^«p: 



(29) 



with 



E^[l- (2(7£p/3)2]i/2 , E' = [{2q£p/if ~ l]^'^ . (30) 

We stress that all these analytical results Eqs. ([TT|) - ([T3l) . 
(HH), ((ni), (imi-dSni) are only apphcable if excluded vol- 
ume effects are negligible. When we consider very long 
semiflexible chains, such that Up ~ L/ ip > n*{ip), we ex- 
pect that the Gaussian results (i?^) = 2£pL — 2£pnp and 
= {l/3)£pL = (1/3)^2 „p hold roughly up to n;{£p), 
and there a smooth crossover to the excluded volume 
power laws, Eq. occurs. We first note that hence 

IT'pi^p, 



corresponds to a crossover radius R* of the chains 



as weh, R*^ 
unity, we conclude 



2ipL* = 2£pn*. Omitting factors of order 



R* 



(31) 



where in the last step Eq. (|20|) was used. For Up > n*{£p) 
we hence expect, invoking the fact that the crossover in 
the linear dimensions for Up = n* should be smooth. 



{R')=R*'{np/n;r 

oitjB'iDllpf-nl- 



,6/5 



(32) 



where in the last step the Flory estimate « 3/5 was 
used (recall that in Eq. ((3T|) the exponent C defined above 
has also been put to its Flory value, C = 2). In terms of 
iV and £p, Eq. ^ becomes (i?^) oc ll'^ {ipDf/^ N^'^ . 
In terms of the constant t^' defined in Eq. (ITil) , we would 

have^^cx^^^/^(^pL')2/5. 

The consequences for the scattering function S{q) are 
now clear, since the gyration radius shows the same scal- 
ing behavior as (i?^), apart from prefactors of order unity. 
Hence we have 



^)^il'\ipDf/'>N^ro 



and only for l/y(R'^) < q < 1/R* we can expect to see 
the nontrivial power law 



S{q) « g-i/-^ 



(34) 



(33) 



while at q* defined from q*R* = 1 we have a smooth 
crossover to the standard Debye law, S{q) oc g"^. Near 
qip = 1 then a smooth crossover to the rod-like scatter- 
ing law S{q) oc q^^ occurs. So the three power laws for 
the radii as a function of chain length (y/ (i?^) cx in 
the rod regime, esc TV^/^ in the regime of Gaussian coils, 
and cx iV" in the regime of swollen coils) find their coun- 
terpart in the scattering function, if TV is large enough. 
The schematic Fig. [5] illustrates these crossover behav- 
iors. The three regimes of the (i?^) versus N (or Up, 
respectively) curve, namely rods, Gaussian coils, and 
swollen coils (Fig. [S^) appear in the S{q) vs. q curve 
(or qLS{s) vs. gL-curve, in the Kratky representation) 
in inverse order: the rods occur for large g, then occurs a 
first crossover to Gaussian coils, and a second crossover 
to swollen coils. Of course, if the chains are very stiff but 
not extremely long, it may be that the regime Up > n* 
is not reached: then in part (a) the swollen coil regime is 
absent, and in part (b) as well: then the K-P model can 
describe S{q) fully, including the regime of the maximum 
of the Kratky plot. Since the crossovers are smooth, it 
may be difficult to identify the different power laws in 
Fig. [5J5 in practice, however. 

We also note that the different regimes are also only 
well separated if both £p is very large (in comparison 
to £b) and also £p/D needs to be very large. If £p is 
very large^but D also (as in the case of bottle-brush 
polymers [11 [13 [zO, IZll) then the regime of Gaussian 
coils disappears from both Fig. [S^ and [SId, and the K-P 
model loses its applicability. 



III. MONTE CARLO SIMULATION METHODS 
AND MODELS 

In the present work, we focus on lattice models exclu- 
sively, because for them particularly efficient simulation 
methods exist; pertinent work on coarse-grained off lat- 
tice models of bottle-brush polymers studied in Molecular 
Dynamics methods for variable solvent quality (98| will 
be mentioned in the conclusions section. 

The archetypical lattice model of a polymer is the self- 
avoiding walk on the simple cubic lattice [o^. Each ef- 
fective monomer takes a single lattice site, the length 
of an effective bond is the lattice spacing, so adjacent 
monomers along the chain are nearest neighbors on the 
lattice. Double occupancy of lattice sites being forbid- 
den, excluded volume interactions under very good sol- 
vent conditions are modelled. 

The properties of this basic model are very well es- 
tablished 100]. Solvent quality can be included as a 
variable into this model implicitly, by allowing for an 
(attractive) energy e that is won if two monomers (that 
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FIG. 5. (a) Schematic plot of the normahzed mean square radius (R^) /{2£pL) versus Up = L/£p (apart from a factor of 2 
this is the number of Kuhn segments), on log- log scales. The Kratky-Porod (K-P) model describes the crossover from rods 
{{R^) = L^) to Gaussian coils ((-R^) — 2£pL). At rip = {£p/D)^, according to the Flory theory a crossover to swollen coils 
occurs, where (R^) oc with v — 3/5 (according to the Flory theory), (b) Schematic Kratky plot of the structure factor of a 
semiflexible polymer, qLS{q) plotted vs. qL, on log-log scales. Four regimes occur: in the Guinier-regime, S{q) ~ 1 — q'^ (R^) /3; 
it ends at the maximum of the Kratky plot, which occurs roughly at qmax-\/ (-R^) ~ 1 (constants of order unity being ignored 
throughout). For very large L then a regime of swollen coils with S{q) oc q"^^'^ is observed, until near qR* ~ 1 a crossover to 
Gaussian coil behavior occurs {R* ^ £p/D). In the Gaussian coil regime S{q) cx q~^, until at q£p of order unity the crossover 
to the rod- like regime occurs {qLS{q) = vr). Oidy the latter two regimes are captured by the Kratky-Porod model. 



which would mean an estimation of 



are not nearest neighbors along the chemical sequence 
of the chain) are nearest neighbors on the lattice. One 
then finds that the Theta point, at which (apart from 
logarithmic corrections 0, lei) the mean square radius 

(R^) scales like a Gaussian chain, (i?^) cx N, occurs with TVpair the number of nonbonded nearest neighbor 



ZN{q,qb) 



config, 



CAr,Ar,.„„iVp.,'Z^^"^g'^-'' , (36) 



for q EE exp(e/fcBT) = qe {- 1.3087) [H. On the 
other hand, if one introduces an energy cost Sb when- 
ever the walk makes a turn by ±90° (of course, reversals 
by 180° are forbidden, because of the excluded volume 
constraint), one can vary the local intrinsic stiffness of 
the chain (cf. Fig. HJ:, which illustrates this model for 
d — 2 dimensions). For qb — exp(— Eb/fesT) = 1 one 
recovers the standard SAW, while the limit — > oo cor- 
responds to rigid straight rods. Following up on our pre- 
vious work [nl, m, [to, [zH 113, we shall focus on this 
model in the present paper, a pplying t he pruned-enriched 
Rosenbluth method (PERM) [i3,[lOl|- PERM is a biased 
chain growth algorithm with resampling and allows to get 
accurate data up to iV = 50000 for this model [t^, 
PERM yields a direct estimate of the partition function 
of a self-avoiding walk with N steps and A^bend 90°-bends 



ZNiqb) - ^ (^'N,N^,,„aqb 

config. 



(35) 



where CN,Ni,^ad the number of configurations of SAW's 
with N bonds and a number A^bond of ±90° turns. It 
would be interesting to extend the approach from ather- 
mal semiflexible chains {q = exp{e/kBT) = 1) to semi- 
flexible chains in solvents of variable quality (q > 1), 



pairs of monomers in the considered configuration. How- 
ever, we are not aware of any study of the full problem, 

Eq. dsni), yet. 

Sampling suitable data on the monomer coordinates of 
the configurations that contribute to the partition func- 
tion Eq. (I35|) . one can obtain reasonably accurate esti- 
mates of the radii and of S{q), as defined in Sec. 2. 

The second model that is studied here is the bond fluc- 
tuation model of bottle-b rush polymers. In the bond fluc- 
tuation model (io3 - [To3 |. each effective monomer blocks 



all eight corners of the elementary cube of the sim- 
ple cubic lattice from further occupancy. Two succes- 
sive monomers along a chain are connected by a bond 
vector £b, chosen from the set {(±2,0,0), (±2, ±1,0), 
(±2,±1,±1), (±2,±2,±1), (±3,0,0), (±3,±1,0)}, in- 
eluding also all permutations. Originally configurations 
were relaxed by an algorithm where a monomer of the 
chain is chosen at random, and one also randomly chooses 
one of the six directions (±a;, ±y, or ±2:), respectively, 
and attempts to move the monomers by one lattice unit 
in the chosen direction. Of course, the move is accepted 
only if it does not violate excluded volume or bond length 
constraints. This move is called the "L6" move. Recently 
Wittmer et al. @ provided evidence that a much faster 
algorithm results if one allows monomers to move to one 
of the 26 nearest and next nearest neighbor sites sur- 
rounding a monomer. With this "L26" move bonds can 
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cross one another, and while such moves do not corre- 
spond to a real dynamics of macromolecules, it leads to 
a much faster expl orati on of phase space and hence a 
faster equilibration |l05j |. 

This model for linear polymers is generalized to the 
bottle-brush architecture by adding side chains at regu- 
lar spacings l/u (which must be integer, e.g. for cr = 1/2 
a side chain is attached to every second monomer of the 
backbone; the densest packing that is studied here is 
cr = 1). The side chains have chain length Ng, and are 
described by the bond fluctuation model as well. Further- 
more, one more monomer is added to each chain end, to 
clearly identify the latter. The number Nh of monomers 
that constitute the backbone then is related to the num- 
ber of side chains ric via Ni, = {ric — 1)/<j + i and the 
total number of monomers of the bottle-brush polymer 
is iVtot = Nb + TLcNs- For the sake of computationa l efln - 
ciency, the L26 move is combined with Pivot moves |lOOl | . 
We refer to [13, llOSf for implementation details. 

As an example for well-equilibrated bottle-brush poly- 
mers as studied in [ll|j IS IzO, [Zl| and in the present 
paper, Fig. [5^ shows selected snapshot pictures for side 
chain length Ns — 24 and various backbone chain lengths 
Nb. According to the visual impression, it seems rather 
natural to describe these bottle-brush polymers by the 
worm-like chain model, but as we shall see below, this 
conclusion would be totally misleading. Experimental- 
ists often are led to a similar conclusion from microscope 
images of semiflexible polymer s (e.g. DNA) adsorbed at 
a substrate (see e.g. jl06l - ll08j |). However, such a con- 
clusion is misleading for several reasons: (i) depending 
on the speed of adsorption of the polymer on the sub- 
strate, the conformation of the adsorbed polymer may be 
a frozen "projection" of the three-dimensional coil, which 
did not have enough time to relax to the two-dimensional 
equilibrium, (ii) In d = 2 dimensions, excluded volume 
forces render the Kratky-Porod (K-P) model of worm- 
like chains inapplicable |62l . l66j , one encounters a direct 
crossover from the rod regime to two-dimensional self- 
avoiding walk behavior (cf. Eq. (jl9p ) when the contour 
length L exceeds the persistence length £p. One also 
should note that the persistence length of a polymer in 
d = 2 dimensions is not at all identical to the persistence 
length of the same polymer in d = 3 d i men sions [63 , [66| . 
The experimental work (see, e.g., 106l - ll08l |) seems to be 
unaware of these problems and the resulting conclusions 
from this work need to be considered with care. It is 
also interesting to note that the snapshot pictures of the 
semiflexible SAW model (Fig. [S]:) do not yield an imme- 
diate visual impression that the chains can be described 
by the K-P model, because of the 90° kinks; however, as 
we shall see, despite this difference of the local structure 
the statistical properties on the mesoscopic length scales 
are well described by the K-P model, for < 10^^ and 
Up less than n*{£p). Thus, we argue that on the basis of 
the inspection of AFM images of semiflexible polymers 
one should be very careful on drawing conclusions which 
model is appropriate to describe these polymers. 
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Nb = 387 
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Nb = 643 



(b) 
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Nb = 195 



Nb = 643 
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= 0,05, N = 25600 



= 0.005, N = 25600 



FIG. 6. (a) Snapshot pictures of bottle-brush polymers as 
described by the bond-fluctuation model, for side chain length 
Ns = 24, and backbone chain length Nt = 99, 195, 387, 
and 643. (b) Same as (a), but displaying the backbone of 
these bottle-brush polymers (c) Snapshot pictures of the SAW 
model with = 25600 and two choices of qt, qt = 0.05 and 
0.005. 



IV. SIMULATION RESULTS 



A. Mean square gyration radii and their analysis 



We start with a description of our results for the 
mean square gyration radius (Rg) of the semiflexible 
SAW model (similar data for the mean square end-to- 
end distance (i?^) of this model have already been pre- 
sented elsewhere [zO, [zH), Fig. [71 We clearly see that 
there are three regimes (Fig. [7^): in the chosen normal- 
ization where we divide out the asymptotic power law 
(i?^) cx N"^", we flrst have a regime where SiRD/ibN^" 
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FIG. 7. (a) Log-Log plot of the relaxed mean square gyra- 
tion radius 3{Rg) /ibN'^'^ versus N, for chain lengths A'^ up 
to A'^ = 50000, and many values of the stiffness parameter 
qb, as indicated. The straight line with slope 2 — 2i/ shows 
the slope reached for small A'^ in the rod-like regime, where 
(iJg) = £lN^ /12; the straight dotted line with slope 1-2:/ 
(for intermediate values of A'^) indicates the behavior expected 
for Gaussian chains, {Rg) = £b£pN/3. (b) Log-log plot of 
3{R^) /£pL versus Up = N£b/£p, using the same data as in (a) 
to test the K-P model (full curve). The different choices of qb 
are shown by different symbols, as indicated. The slope that 
one expects for all qb for A' — s- co , (21/ — 1), is indicated by a 
broken straight line. 



increases with N'^~'^^ . For small qi, this is interpreted as 
a rod-like regime; for qi, > 0.2 the chains are still too 
flexible, however, so a strictly rod-like behavior cannot 
yet be seen. Then a maximum occm'S, and the ratio 
3{Rg) /£bN'^'^ decreases, before it settles down, after a 
second smooth crossover at a horizontal plateau (which 

according to Eq. (|14p defines the value i^^). While this 
plateau for qi, > 0.05 is (presumably) actually reached 
for N — 50000, the data also indicate that for qi, < 0.03 
even chains of length iV=50000 are at least an order of 
magnitude too short to allow a direct convincing esti- 
mation of the amplitude value £p " . On the other hand, 
even for qi, = 0.005 (where we estimate from Eq. H) that 
the persistence length £p is as large as £p ss 52 [B^) the 
slope of the data in the intermediate regime has not fully 
reached the theoretical value 1 — 21/, the slope of the data 
in Fig. [7^ is still affected by crossover effects: the gradual 
crossover away from the Gaussian plateau towards the ex- 
cluded volume-dominated behavior already starts when 
the gradual crossover from the rod-like regime to the 
Gaussian regime ends. Thus, even stiffer chains would 
be required to have a fully developed Gaussian behavior 
of the gyration radius. Fig. [7]d now attempts a scaling 
plot, where the persistence length estimates extracted 
from Eq. ^ were used to rescale (Rg) in the K-P model 
representation (cf. Eq. (fT7)) ). It is evident that the rod- 
like regime and the onset of the crossover towards the 
K-P plateau are very well described by Eq. (flT)) . For 
qb > 0.2, of course, there is basically a direct crossover 
from the rod-like regime to the excluded volume domi- 
nated regime, but even then it is evident that the curves 
do not superimpose on a master curve, as they do in 
d = 2 dimensions [g^, , but rather splay out systemat- 
ically, and the smaller becomes (and hence the larger 
£p becomes) the more the data still are slightly above the 
K-P plateau. 

Using the estimates for n*{qf)) extracted from the anal- 
ysis of (i?^) for this model in our previous work (toI. l7l|. 
the data for Up ^ n* do collapse on a simple straight 
line on the log-log plot, however (Fig. [T]:). For Up near 
n* the curves splay out, the master curve describing this 
second crossover from the K-P plateau to the excluded 
volume power law emerges as an envelope of the curves 
for individual values of qt (which fall increasingly below 
the master curve in the crossover region the larger qi, 
is). Of course, n* cx {tp/tbY cannot produce a scaling of 
the crossover towards the rod-behavior, there the curves 
must splay out, irrespective of how small is, but the 
deviation of the data from the horizontal K-P plateau 
moves more and more to the left of the plot the smaller 
qh becomes. 

Recalling that for the semiflexible SAW model the ef- 
fective chain thickness D simply is I? = ^f, = 1, the Flory 
theory, Eq. ([20)) . simply predicts n* oc ip] ( — 2, while 
the rod to Gaussian coil behavior occurs around Hp = 1, 
of course. Qualitatively, our data are in good agreement 
with these predictions, but not quantitatively: This is 
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FIG. 8. Log-log plot of n* and Ip" /Hb, as indicated in the 
figure, versus Ip/ib, using data for qt = 0.2. to 0.005 (left to 
right). Dotted lines indicate the observed exponents while full 
straight lines show the Flory predictions for the exponents. 



illustrated in Fig. El where n* and Ip^ are plotted in 
log- log form versus ip/tb- It is seen that instead of the 
theoretical value ^ = 2 an exponent ^ = 1.5 is observed. 
Now it is clear that Flory arguments imply also = 3/5 
instead oi v 0.588 '7], but this small difference cannot 
account for the large discrepancy encountered here. It 
would be desirable to study much larger values of N to 
confirm whether this discrepancy is a real effect (or our 
estimation of the crossover master curves in Fig. [7h for 
(and for {R^) in 71]) are systematically off). Thus, 
more work is still needed to fully clarify the situation. 

It turns out that the behavior of our model for the 
bottle-brush polymers (which can describe actual scatter- 
ing data for bottle-brush polymers very well, as demon- 
strated by Hsu et al. [13]) is much simpler: a plot of 
the mean square gyration radius (i?^ of the backbone 
versus backbone chain length N^, for different side chain 
lengths Ns, Fig. [5^, normalized by A^^" reveals a mono- 
tonic increase towards a plateau, there is not the slight- 
est indication of a regime where the ratio (i?^ /N"^" 
decreases, unlike the behavior of the semiflexible SAW 
(Fig. [7^). Thus, there is no evidence whatsoever for a 
Gaussian K-P plateau for this model. But the increase 
of the plateau value Ip " with increasing side chain length 
Ns does indicate that the chain considerably stiffens, as 
Ns increases. However, this stiffening goes along with an 
increase in the effective chain thickness D. The latter 
can be estimated from the radial density profile (Fig. [TU)) 
by identifying the diameter D of the bottle-brush as 
D = 2rJ{Ns) = 2/(i?2j. Hsu et al. suggested 
to coarse-grain the bottle-brush along the backbone, di- 
viding it into "blobs" along the chemical sequence of the 
backbone. The chemical distance Sbiob along the back- 
bone between its exit and entry points into a blob is 
found from a simple construction which assumes that 
the blobs are essentially spherical, so the geometrical dis- 
tance Ar(s) between exit and entry points of the back- 
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FIG. 9. (a) Plot of the normalized mean square gyration ra- 
dius of the backbone of the model for bottle-brush polymers, 
3{R^ fj)/^bN'^'^, versus the backbone chain length Ni,, for four 
difi'erent side chain lengths Na, Ns = 6, 12, 18, and 24, as 
indicated. Here it ~ 2.7 is the average bond length of the 
single bond in the bond fluctuation model under good solvent 
conditions. The horizontal plateaus allow to extract the esti- 
mates for Ip" defined in Eq. H14|) . (b) Same data as in (a), 
but ordinate is rescaled with £f'' so for Nt oo all data con- 
verge to one, and abscissa is rescaled by the effective blob-size 
Sbiob (see text), on a log-log plot. Straight line indicates the 
rigid-rod behavior, with slope 2 — 2i/ in this representation. 



bone should be equal to D. Recording Ar(s) for arbitrary 
s, Fig.lTUb. using the equation Ar'(sbiob) = D allows us to 
simply read off the numbers Sbiob for the choices of Ns , as 
illustrated in Fig. (TU] The success of the rescaling shown 
in Fig. [Hb shows that the persistence length ip of bottle- 
brushes simply is proportional to D. We also recognize 
that the asymptotic SAW-like behavior (where the hori- 
zontal plateau in Fig.lHlis reached) only occurs for about 
Nb/shioh ~ 60. Comparing this behavior to Fig. [71d, we 
see that there the power law (for qb = 0.4) or K-P plateau 
(for small qb) is reached for Up — N£b/£p « 40. Roughly, 
these successive blobs then are equivalent to one per- 
sistence length. This comparison suggests that for the 
bottle-brushes we should identify |(iVb/sbiob) with Up, 
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FIG. 10. (a) End-to- End distance Ar(s) of subchains con- 
taining s successive backbone monomers for Ns — 6, 12, 
18, and 24. The horizontal sohd hnes and the numbers 
shown on the ordinate indicate the choices Ar-(s) — 2Rc3{Ns), 
from which the corresponding values of Sbiob can be read 
off (vertical straight lines), namely, Sbiob = 6, 10, 12, and 
14 for Ns = 6, 12, 18 and 24, respectively. Rcs{Ns) is 
the cross-sectional radius, which is extracted from the radial 
monomer density profiles p(r) as shown in (b). (b) Radial 
monomer density profiles p(r) in planes locally perpendicu- 
lar to the backbone of bottle-brush polymers with backbone 
length Nb = 1027 and plotted versus radial distance r for side 
chain lengths Ns = 6, 12, 18 and 24, as indicated. The cross- 

oo 

sectional radius then follows as (i?cs) = 27r J rdrp(r)r^ with 



the density profile being normalized as 2-n J rdrpir) — Ns. 



The values of Res are pointed out by arrows. 



i.e. the number of monomers along the backbone corre- 
sponding to one persistence length is |sbiob = 9, 15, 18 
and 21 for Ng = 6, 12, 18 and 24, respectively. Noting 
that the average bond length £b in the bond fluctuation 
model is £i, — 2.7, we would obtain persistence lengths 
ipiNs) = K 17, 27, 34 and 41 for N, = 6, 12, 18 
and 24, respectively. The result that ip{Ns) is of the 
same order as D irrespective of the side chain length 
agrees with early theoretical predictions [l^ 13 but is 



at variance with the result of Fredrickson [20j who pre- 
dicted a much faster increase of ip with Ng. However, 
Feuz et al. [13] pointed out that the result of Fredrick- 
son [2^ can only be expected to hold for extremely long 
side chains, such as Ns — 1000. Such long side chains 
are neither relevant for simulations nor for experiment, 
however. We stress that the range of Ng accessible to 
simulations (Figs. [SI fTU)) nicely corresponds to the range 
of studied experimentally [H, [s^ - lssj . 

The mapping performed in Fig. [S] means that we have 
coarse-grained the bottle-brush polymers (Fig. into 
an effective bead-spring model (Fig. |3Jl). If this map- 
ping is taken literally, it can also be used to obtain the 
resulting coarse-grained contour length Lcc (Fig. |31) as 



Lr 



2R,,{Ns)Ni,/s,,iohiNs). 



(37) 



Instead of the "chemical" contour length L — Ni,£i, w 
2773 a reduced length is found, namely Lcc ~ 1989, 1824, 
1963 and 2016, for Ns = 6, 12, 18 and 24, respectively. 
This means that the coarse-grained contour length Lcc is 
about 30% smaller than the "chemical" contour length 
in this model. 



B. Analysis of the structure factor 

We now turn to the structure factor of the semiflexi- 
ble SAW model presenting Kratky plots for two contour 
lengths, L = 400 and L = 25600, in Fig.HH As expected, 
cf. Fig. O one first has a linear increase with qL, then a 
round maximum followed by a decrease which then grad- 
ually settles down at a horizontal plateau, that again is 
compatible with the theoretical prediction, tt. 

While for the short chain length {L = 400) the agree- 
ment with the theoretical prediction (due to Kholo- 
denko 89], Eqs. (I2Hl)-(ISni), which were found 97] to be 
numerically almost indistinguishable from the exact re- 
sult provided by Stepanow [H, is almost perfect, 
for the very long chains {L = 25600) we note system- 
atic deviations between Kholodenko's prediction [80| and 
the data for relatively large qb near the maximum of the 
Kratky plot. This must be expected, since the input in 
the Kholodenko formula is just the persistence length ip 
{which we have extracted from Eq. ^} and implicit in 
the theory is the Gaussian prediction for (Rg), namely 
(Rl) = Lip/3 {Eq. dni)}. As seen in Fig. El for qt = 0.2 
already rather clear deviations from this result occur for 
N = 25600, while for small N such as iV = 400 such 
deviations still are rather small. In contrast, in the De- 
bye formula the correct (as observed) value of (Rg) was 
used as an input: then deviations from the Debye func- 
tion are only seen near the region where the crossover to 
qLS{q) — TT starts to set in at large q (the Debye function 
does not describe this crossover at all). Since the shape 
of the Kholodenko function always is rather similar to 
the actual function, it is obvious that one always can fit 
the data to the Kholodenko function, if ip is not known: 
however, the resulting fitted persistence length will be 
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FIG. 11. Kratky plot of the structure factor of the semiflexible 
self-avoiding walk model, log-log plot of qLS{q) versus qL, for 
two chain lengths, namely L = 400 (a) and L = 25600 (b). 
(c) Persistence length £p plotted vs. qt for 0.005 < qt < 0.4. 
In (a)(b) several choices of qt are included, namely qb = 0.2, 
0.05, 0.02, 0.01, 0.005, as indicated. The scattering functions 
of a rigid rod and the Debye function are included, as well as 
the prediction of Kholodenko {Eqs. ([281) -dSO])}. The predicted 
large q-limit of tt is indicated. In the Debye function the 
observed value of (-Rj) was used as an input, while for the 
Kholodenko formula the persistence length (estimated from 
Eq. Q) was used as an input. In (c) the persistence length is 
taken as the best fitting parameter such that the prediction of 
Kholodenko formulas describes the correct maximum in the 
Kratky plot for our simulation data of chain lengths L — 400, 
and 25600. The estimates using Eqs. ([l}, ^ are also shown 
in (c) for comparison. 
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FIG. 12. Kratky plot for qt = 1.0, 0.2, and 0.005 plotted in 
a rescaled form, qS{q)/ {qS{q))nis.x vs. q/gmax, restricting the 
ordinate range to the decade from 1 to 0.1 and the abscissa to 
the range from 1 to 10. The theoretical power laws (slope= 
1 — l/u « 2/3 and slope — —1, respectively, as predicted in 
Fig. 5) are included. 



systematically too large, if excluded volume effects are 
present as shown in Fig. [TTh . 

To elucidate the significance of excluded volume on the 
structure factor further, we show a magnification of the 
region near the maximum for qt — I, 0.2, and 0.005 in 
Fig. [121 It is seen that the identification of the two power 
laws suggested for the decay ofqS{q) in the region beyond 
the maximum of the Kratky plot is rather subtle. In par- 
ticular, for rather stiff chains the crossover to the rod-like 
scattering sets in rather early, so for the clear identifica- 
tion of power laws the available range of q simply is not 
large enough. This very gradual crossover between the 
three different regimes (rods to Gaussian coils to coils 
swollen by the excluded volume interaction) complicates 
the data analysis, if only a restricted range of parameters 
(such as the chain length N and the wavenumber q) can 
be investigated. 

The smoothness of the crossover also becomes evident 
when one studies the dependence of the position of the 
peak in the Kratky-plot (and its height) on the persis- 
tence length (Fig. [T3|) . Typically, the data fall neither in 
the regime where strict Gaussian behavior occurs, nor in 
the regime where excluded volume scaling is fully devel- 
oped. 

Despite all these difficulties due to the gradual 
crossovers, the semifiexible SAW nevertheless is a rela- 
tively simple case, since one knows that here D = £b{= 
1), and £p can be varied over a wide range by variation 
of qb, keeping all other parameters constant, and more- 
over £p can be estimated precisely from the initial decay 
of the bond vector autocorrelation function (or, equiv- 
alently, from Eq. (jj])). For the second model studied 
here, bottle-brush polymers under good solvent condi- 
tions, we have seen that varying the side chain length Ng 
we change D and ip together, and also the coarse-grained 
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estimated using Eq. Q) using data for qb = 0.2 to 0.005. 



Straight lines indicated the exponents 

-1/2 



oc L 



-1/5 



and 



gmax oc £p that one expects according to the Flory treat- 
ment in the excluded volume region and Gaussian region, re- 
spectively. All data were taken for N=50000. 



contour length Lcc is significantly smaller than the chem- 
ical contour length Nhih, and it is nontrivial to estimate 
Lcc accurately. 

In previous work [ill [iTj we have already considered 
the decomposition of the total scattering function of 
bottle-brush polymers into the scattering from the back- 
bone and from the side chains. This analysis which has 
the advantage that it provides a direct link to correspond- 
ing experiments [H, [s^] will not be addressed here, but 
we rather focus on the scattering function of the back- 
bone only. Fig. [U shows Kratky plots for relatively 
short backbone chain lengths {Nf, = 131 and 259, respec- 
tively). One recognizes that for short side chain lengths 
{Ns = 6,12) qLSb{q) does not settle down to a well- 
defined "Holtzer Plateau", at least not within the avail- 
able window of wavenumbers. Clearly, also the range over 
which qLSb{q) decays from the maximum to the horizon- 
tal part that appears for Ng = 24 and 48 is rather small, 
and does not warrant any analysis in terms of the power 
laws suggested in Fig. [S]d. This mismatch between the 
actual plateau values (for Ng = 24 and 48), which are 
close to 4, and the theoretical value n can be attributed 
to reduction of Lcc in comparison to Nhit since the ac- 
tual orientations of the backbone vectors are not strictly 
aligned with the coarse-grained backbone (Fig. 2]), as is 
also evident from the fact that {Si ■ ai+i) / {af) (Fig. [2)d) 
is already reduced to about 0.7, but the further decrease 
of {di ■ (ii+s) / {of) is rather slow, due to the side chain 
induced stiffening of the backbone on mesoscopic scales. 
The ratio at about 4/7r is compatible with the reduction 
of Lcc by about 30% relative to L noted previously, so 
gratifyingly our analysis is internally consistent. 

Fig. [12] shows plots of S{q) vs. q for fixed side chain 
length Ns=2A but different backbone chain lengths. In 
this plot, an attempt is made to locate an onset wavenum- 
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FIG. 14. Kratky plot of qLSb{q) versus qL for bottle-brush 
polymers with A*';, — 131 (a) and Nb — 259 (b), where 
"L = Nbib" is the "chemical" contour length and Sb{q) is 
the scattering function of the backbone only. Four side chain 
lengths Na = 6, 12, 24 and 48 are included, as indicated. 
The horizontal straight line shows the Holtzer plateau (= tt) 
if the coarse-grained contour length could be identified with 
the chemical contour length. 



ber q* for the Holtzer plateau, in terms of a fit of two in- 
tersection straight lines. Of course, the data are smooth 
and the onset of the Holtzer plateau does not occur 
sharply but rather gradual; thus q* can be estimated 
only with considerable error (for large Nb we estimate 
q* w 0.06 ± 0.01, while for Nb = 131 the estimate rather 
is q* sa 0.075 ±0.020 II]). Now the question is, how can 
one relate q* explicitly to the persistence length? Should 
one take ip — 27r/g*, or ip = l/<?*? Lecommandoux et 
al. [l^l who were the first to try such a method suggested 
the relation ip « 3.5 /q*, but we see little theoretical sup- 
port for this choice either. 

It would be advantageous if one could rely on the des 
Cloizeaux relation, Eq. (^5]) . which suggests to plot qS{q) 
vs. 1/q for q£p ^ 1: one should find a straight line, the 
intercept at the ordinate should yield tt, the slope of the 
straight hue should yield 2/{3£p). 

However, when one tests this method for the semiflex- 
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FIG. 15. Log-log plot of qSb{q) vs. q for bottle-brush poly- 
mers with side chain length Ns = 24 and various backbone 
chain lengths Nb from Nt = 131 to Nb = 1027, as indicated. 
The point of intersection between two Broken straight lines 
illustrates the estimation of q*{Nb), the wavenumber where 
the onset at the Holtzer plateau occurs. 



ible SAW, one finds that the data that can be fitted to a 
straight line are at q£p « 1 rather than at qip ^ 1, and 
the slope of the straight fine disagrees with the predic- 
tion (Fig. 116b.). Thus, it is not really a big surprise that 
this does not work well for our bottle-brush model either 
(Fig.IISlD). 

An interesting alternative of data analysis is, however, 
a fit of the Kholodenko formulas {Eqs. (Il8l)-(l3Dl)} to the 
structure factor, using both Lcc and £p as individual ad- 
justable parameters for each value of Nb (Fig. fT7|). First 
of all, one sees that the Kholodenko structure factor pro- 
vides a good fit in all cases, and the result for the coarse- 
grained contour length, Lcc — 2.03Ni„ even is physically 
very reasonable: we have obtained that there is a 30% 
reduction of the Lcc in comparison to the "chemical" con- 
tour length L = Ni,li, — 2.7N\j in the previous subsection. 

However, the problem of this fit is the unphysical be- 
havior of the persistence length £p: since we know that 
the Kholodenko [s^ approach involves necessarily the 
Gaussian result (i?^) = \tpLcc but we know that for 



our model (i?^) cx N'^" and Lcc ~ 2.03A^f,, the only way 
to reconcile these results is a persistence length scaling 
as ip cx N^'^~^, and this is what we see in Fig.lTTb. Thus, 
despite the seemingly good fit (Fig.fTTk) and good results 
for Lcc fFig. 117b). the results for the persistence lengths 
are completely unreliable! 

In order to apply this approach, one must make sure 
that one works with data in the Gaussian regime, and 
this is not at all the case for bottle-brush polymers under 
good solvent conditions. 
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FIG. 16. (a) Plot of qLS{q) vs. for the semiflexi- 

ble SAW and L — 25600, including 4 choices of qb, as indi- 
cated, (b) Same as (a) for bottle-brush polymers with fixed 
side chain length Ns — 24 and varying backbone length as 
indicated. Broken straight line is the des Cloizeaux 80] pre- 
diction, Eq. (I26p . full straight line an empirical fit to the data. 



V. CONCLUSIONS 

In this paper, we have focused on the behavior of sin- 
gle semiflexible polymers under very good solvent con- 
ditions, considering how the chain stiffness affects poly- 
mer properties such as the mean square gyration radius, 
the structure factor, etc. Our analysis focused on the 
question how the variation of chain stiffness affects these 
properties, and hence one can infer from these properties 
a characterization of the "intrinsic stiffness" of the poly- 
mer chain in terms of the so-called "persistence length" . 

We have contrasted two models, the self-avoiding walk 
on the simple cubic lattice where a bending energy 
Eb causes pronounced stiffening of the polymer when 
Eb ^ ksT, and a lattice model for bottle-brush poly- 
mers, where backbone stiffening is caused by increasing 
the length of side chains. These two models constitute 
two quite distinct limiting cases: in the SAW model, in- 
crease of Eb/ksT causes stiffening without any effect on 
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FIG. 17. Plot of qLccSb{q) versus Lcc? for bottle-brush poly- 
mer with side chain length A'^s = 24 and several values of Ni,, 
from Nb = 131 to 1027, as indicated. For each choice of Nt 
both Lcc and £p were individually adjusted, (b) Plot of Lcc 
vs. A^i,, showing that Lcc ~ 2.03A^i, rather than being equal to 
the "chemical" contour length L = ItNt « 2.7Nb. (c) Plot of 
the persistence length £p (from the fit in (a)) versus A*';,. In 
the log-log plot (inset), the theoretical power law £p oc N^'^~^ 



is shown by the straight line. 



the local thickness of the chain, which strictly remains the 
lattice spacing. For the bond fluctuation model of poly- 
mer brushes, however, we have found that backbone stiff- 
ening is caused by the thickness of the (coarse-grained) 
worm-like chain, the persistence length increases propor- 
tional to the cross-sectional diameter of the bottle-brush. 

Since snapshot pictures (Fig. ^ suggest that the 
bottle-brush polymers (or their backbones, respectively) 
resemble worm-like chains (and the same conclusion is 
often drawn from AFM pictures or electron micrographs 
of actual polymers) , the use of the Kratky-Porod worm- 
like chain model has become very popular. However, 
we demonstrate here that for bottle-brush polymers this 
model yields very misleading results: since the mean 
square gyration radii of bottle-brushes are found to scale 
with their contour length L as (Rg) oc L^'^, the Kratky- 
Porod (K-P) result (Rg) — £pL/3 invariably causes a 
spurious contour length dependence of the persistence 
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length when fit to the data, namely ip{L) oc L 
as L oo. Although the fits of the K-P model look al- 
most perfect (Fig, 117b ) and numbers for L resulting for 
the contour length from the fit are rather reasonable, the 
result for "the" persistence length simply is meaningless! 

Already in our earlier papers we have shown that sim- 
ilar ambiguous results for the persistence length are got- 
ten when orientational correlations along the chain back- 
bone are analyzed, or the projection of bond vectors on 
the end-to-end distance are studied (although the result- 
ing numbers for £p{Nb) seem to be somewhat smaller 
than those shown in Fig. [T7b). The large g-behavior of 
the structure factor S{q) yields a qualitatively more rea- 
sonable behavior, but a unique choice for a well-defined 
persistence length as a measure for intrinsic chain stiff- 
ness does not emerge. All these difficulties in under- 
standing the stiffness of bottle-brush polymers in good 
solvents are intimately linked to the fact that one can 
coarse-grain into some effective self-avoiding walk model 
(Figs, ini nop , and no regime exists where the polymers 
resemble Gaussian chains. Of course, this fact is differ- 
ent if we would consider bottle-brush polymers in Theta- 
solvents (as done by Theodorakis et al. [98]), since then 
oc L and the use of the K-P model is qualitatively 
reasonable. Another interesting possibility to extract 
a persistence length of bottle brushes would be an at- 
tempt to estimate an effective bending modulues. One 
would have to estimate the coarse-grained free energy of 
bent versus non-bent configurations of suitable pieces of 
bottle-brush polymers, which in principle can be deduced 
from sampling suitable angular distribution functions for 
such sub-chains. However, the implementation of such 
an approach is not straightforward and has not been at- 
tempted. 

We have found that the situation in some respects is 
simpler if one considers polymers where the stiffness can 
be enhanced while keeping their thickness constant, as 
modeled by a semiflexible extension of the standard SAW 
model. Then an intermediate Gaussian-like behavior of 
the mean square radii and the structure factor emerges, 
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and this can be understood theoretically (Fig. [5]), at least 
in qualitative terms. While still the asymptotic decay of 
bond vector autocorrelation functions is unsuitable to in- 
fer anything about the intrinsic stiffness (due to the fact 
that the asymptotic decay is not exponential but rather 
described by a power law), in favorable cases the initial 
decay of these autocorrelation functions provided use- 
ful estimates of the persistence length, which then can 
be used as input in the K-P model. While still some 
problems occur to understand for very long chain the 
crossover between the K-P model and the ultimate SAW 
behavior, quantitatively, in qualitative terms the situ- 
ation is understood. We emphasize, however, that all 
these comments only address the three-dimensional case: 
in c? = 2 dimensions, the K-P model does not work at all, 
and one has a direct crossover from rod-like polymers to 



SAW'S. 

It is hoped that our analysis will help experimentalists 
with a proper interpretation of their data on semiflexible 
polymers. 
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